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Low- lying excited states as well as the ground state of the spin-1 anti- 
ferromagnetic Heisenberg chain with a spin-1/2 impurity are investigated 
by means of a variational method and a method of numerical diagonaliza- 
tion. It is shown that 1) the impurity spin brings about massive modes 
in the Haldane gap, 2) when the the impurity-host coupling is sufficiently 
weak, the phenomenological Hamiltonian used by Hagiwara et al. in the 
analysis of ESR experimental results for Ni(C2H8N2)2N02(C104) contain- 
ing a small amount of spin-1/2 Cu^"^ impurities is equivalent to a more 
realistic Hamiltonian, as far as the energies of the low-lying states are 
concerned, 3) the results obtained by the variational method are in semi- 
quantitatively good agreement with those obtained by the numerical diag- 
onalization. 
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Since Haldane's prediction^^ of the difference between integer-spin and 
half-integer-spin antiferromagnetic Heisenberg chains, the spin-1 antiferro- 
magnetic Heisenberg chain has been the subject of a large number of the- 
oretical and experimental studies. One of the recent topics of this subject 
is the edge effect on the chain. Kennedy^-* found that the open chain has a 
fourfold degenerate ground state composed of a singlet and a triplet which 
we call the Kennedy triplet, in contrast to a unique singlet ground state of 
the periodic chain. The fourfold degeneracy of the ground state, which 
was originally found in the so-called AKLT model^^ with open boundary 
conditions, is considered to reflect the hidden x symmetry in the 
open chain. ^) The hidden symmetry is attributed to the spin-1/2 degrees 
of freedom at edges of the chain. 

Recently the present authors and Harada^^ investigated theoretically 
the impurity-bond effect on the ground state and the low-lying excited 
states of the chain to interpolate between the open-chain and periodic- 
chain cases. They showed that the impurity bond brings about a massive 
triplet mode in the Haldane gap and that the triplet state comprises three 
of the four ground states of the open chain. Miyashita and Yamamoto^^ 
performed the Monte Carlo analysis of the open chain to show that the 
magnetic moments localized around the edges for the Kennedy triplet 
decay exponentially with the decay constant which is about 6 in lattice 
spacing. Hagiwara ei al?'^^ performed the ESR experiment on the spin-1 
linear-chain antiferromagnet Ni(C2H8N2)2N02(C104), abbreviated NENP, 
containing a small amount of spin-1/2 Cu^+ impurities and gave for the 
first time experimental evidence for the existence of the spin-1/2 degrees 
of freedom at the host-spin sites neighboring the impurity. They ana- 
lyzed successfully their experimental results by using the phenomenological 
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Hamiltonian given by 

2 

yp^e _ ^ (^J^SqS^ + JjSqS^ + JzSqS^) , (1) 

where Sq is the spin- 1/2 operator of the Cu^"'" impurity; Si and S2 are 
the spin operators which represent the spin-1/2 degrees of freedom at the 
host-spin sites neighboring the impurity; Jx, Jy, and are the effective 
exchange constants. There is, however, no clear explanation for the origin 
of the anisotropy of the effective exchange interaction. 

In this letter, we investigate the spin-1/2 impurity effect on the spin- 
1 antiferromagnetic Heisenberg chain from the theoretical point of view. 
Our final goal is to analyze the above ESR experiment by comparing the 
theoretical results with the experimental ones. As the first step to the goal, 
we discuss the relation between the phenomenological Hamiltonian 7i^^^ 
and the more realistic Hamiltonian given by 

n = no + n' , (2a) 

i=i 1=1 

(J > 0, A > 0) , (2b) 
'H! = J' {sqS^ -\- SqS\ -\- A'sq(S'j) -|- J' {^SqS^ -\- SqS^ -\- X'sqS^^^ 

(A' > 0) , (2c) 

where Sq is the spin-1/2 operator of the impurity spin as is stated above 

— > 

and Si {i=l, 2, • • •, A^) is the spin-1 operator of the host spin. Thus, Ho 
and 7i' represent, respectively, the Hamiltonian for the host-host coupling 
and that for the impurity-host coupling. Using both an analytical method 
and a method of numerical diagonalization, we calculate the energies of the 
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ground state and the low-lying excited states of H.. By comparing these 
energies obtained by the analytical method with those for the Hamiltonian 
we will show that, when | J'|/J is sufficiently small, Ti^^^ is equivalent 
to ?Y, as far as the energies of the low-lying states are concerned. We also 
explore the dependences on J' and on d of the energies of the low-lying 
states, since they yield information on the relation between the results of 
the above ESR experiment and the impurity-host interaction. 

We now calculate the energies of the low-lying states of assuming 
that I J'|/J<Cl, that is, treating Ti' as a small perturbation. In the Haldane 
region, to which we confine ourselves hereafter, the ground state of the 
unperturbed Hamiltonian TIq is fourfold degenerate. The variational 
method discussed in the previous paper, according to which the Haldane 
region is given by 4>c?>2A— 4, leads to the following fourfold ground-state 
wave functions $ and (r = +, 0, — ) expressed in the matrix-product 
form.^'^°) The function $ describes the state with no domain wall and is 
given by 

$= Trace [0102 •• -^A^-i^Tv] , (3) 

sinf6'^ 

= cos{e) Ce (7z + a+ + /3ia-) , (4) 

where a^, and Pi are the spin states at the i-th. site, which correspond, 
respectively, to 5*1 = 1, 0, and —1, and a± [= (ax^icTy) / V^] , ax, cry, and cr^ 
are the Pauli matrices. The parameter 6 is determined from the equation, 

cos(26i) = ^ . (5) 

(t) 

The function ^)^^ describes the states with a domain wall and is expressed 
in terms of 0£ and the wall operator w as 

= Trace [ 0i02 • • • (f>N-i(t>N w] , (6) 
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where w = —a- for r = +, w = Gz for r = 0, and w = (7+ for r = — . It is 
noted that, when A=l and (i = 0, $ represents the singlet state, and 

and represent, respectively, the triplet (Kennedy triplet) states 
with Mo = Xl^i*^! = 1? 0? ^iid ~1- These four wave functions give the 
same energy expectation value £"0,0 and the correlation length ^ as 

^o,o = -(iV-l)j(i-^)(l + J)+iV..^, (7) 



^=-l/ln 



d-\ 



(8) 



4- A 

Performing a perturbation calculation, we restrict the wave functions 
for 7-^0 to $ and ^'^^ . We denote the wave function for sq by Xv^ where 
V = -\- for Sq = +1 and v = — for Sq = — |. Then, the bases of the wave 
functions for the low-lying states of l-i may be represented by the product 
of the former four functions and the latter two functions as 

= ^ X. , ^[."^ = X. ■ (9) 
The matrix elements of TL' in this representation are easily calculated from 

H' ^u = ^ , (10a) 

n' = J' X' sin^ e ^^^^ , (10b) 

n' ^1°) = J' sm{2e) , (10c) 

n' ^^^^ = J' sm{2e) ^1,°) - J' X' sin^ § ^^^^ , (lOd) 

where = + or — , and i' = + when v = — and = — when = +. In 
deriving eq. (lOa)-(lOd), we have neglected the factor cos^ 9. Solving the 
corresponding secular equation, we obtain the energy eigenvalues (mea- 
sured from £"0,0) as 

441) -^'^'W^y (") 



.1 3\ 4-d r / _ 4-2A + d J 



es(±ij)=0. (14) 

Here, we have denoted the energy eigenvalue of the state with M = Sq + 
Yl!i=i^i ^r(-^? S) (r = s, t), where S represents the magnitude of the 
total spin of the corresponding state in the isotropic case of A = A' = 1 
and d = 0. The subscripts s and t show that the eigenvalues are associated 
with the bases and respectively; the wave function for es(i^5 \) 

is given by that for et{-^\, |) by \ and those for et{^\, |) and 
et(i|, \) by linear combinations of '^^^ and 

Before discussing the results of the above analysis in more detail, 
we examine the relation between IHP^^ and Ti. The energy eigenvalues 
eP^«(M, S) (r = s, t) oiUP"^^ is easily calculated to be^) 

1 3 
2' 2 

^f^(±J,^)=-^^Wl + T^ + U, (17) 
.f^(4^)=0, (18) 

where A^fj = 2J^/(J^+Jy), and where M represents the z-component of the 
total spin in the case of Jx = Jy, S represents the magnitude of the total 
spin in the isotropic case of Jx = Jy = X , and the subscript r have the same 



etH±\,h = ^-i{^ll + Tl--l \ , (16) 
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meaning as that of £r{M, S) but for the three-spin- 1/2 system. Compar- 
ing eqs. (15)-(18) with eqs. (11)-(14), we can determine the correspondence 
between the interaction constants in H"^^^ and those in H. The results are 

J. ^ J' y , (19) 

(Ay(4-d) 

2(4 -2X + d) • ^^^^ 

We have thus shown that, H^^^ is equivalent to H when | J'|/J <C 1, as 
far as the energies of the low-lying states are concerned. Equations (19) 
and (20) give a clear explanation for the origin of the anisotropy of the 
exchange interaction in T-^p'^®; as seen from eq. (20), in the case of A = A' = 1 
the uniaxial anisotropy d in 7i produces the anisotropy of the exchange 
interaction in 7i^^^. This result confirms the concept of the spin- 1/2 degrees 
of freedom at edges of the open spin-1 chain. The use of HP^^ for the semi- 
quantitative analysis of the ESR experimental results for the NENPiCu^"^ 
system"^'^) is also justified. 

Let us discuss several qualitative properties of the energies of the low- 
lying states in the case of A = A' = 1, which are deduced from eqs. (11)-(14). 
For convenience, we choose the origin of the energies £^o,o+£^t(=l=|? ^) and 
define Ar(M, S) (r = s, t) as 

Ar(M, 5)=er(M, 5)-£t(±^, ^) . (21) 

In Fig. 1 we show the d-dependence of Ar(M, S)/J'. We see from this 
figure that, when \ J'\/ J is sufficiently small, Ai.(M, S) satisfies the relation 
0<A,(±i, i)<At(±§, i)<At(±i, §) or 0>As(±i, ^)>At(±|, |) > 
At(±|, I) depending upon whether J' > or J' < 0. It should be noted 
that At(d=^, I) is of course equal to At(±|, |) in the isotropic case of 



A = A' = 1 and d = 0, and also that in this case the ratio R defined by 
i? = As(±|, ^)/At(±|, |) is given by |. As has been discussed in ref. 5, 
the magnitude of the Haldane gap, which should be defined as the energy 
difference between the bottom of the energy continuum and the ground 
state, is not affected by the presence of an impurity spin (see Fig. 5 in 
ref. 5). Combining this with the results shown in Fig. 1, we obtain the 
schematic energy versus J' diagram given in Fig. 2. This figure shows that 
the impurity spin brings about the massive modes (the so-called impurity 
states) in the Haldane gap in a certain range of the impurity-host exchange 
constant J'. Figure 2 also suggests that we can determine at least the sign 
J' from the experimental results for energy-level separations. For example, 
when d>0, the Zeeman splitting of the second- lowest energy level due to 
the external magnetic field H will give a key for determining the sign of J', 
because the value of M of the second-lowest energy level at = is either 
±1 or ±1 depending on whether J'>0 or J'<0. 

In order to numerically examine the analytical results obtained above, 
we have performed a numerical diagonalization by the Lanczos method""^^) 
for finite- {N = b, 7, • • •, 15) chains in the isotropic case. In Fig. 3 the 
results for As(±|, ^) and At(±|, §)[ = At(±|, §)] obtained for Ar=15 
are plotted as a function of J'; the dashed line and the dotted line represent 
the former and the latter, respectively. This figure should be compared 
with Fig. 2(b). Numerical results show that At(±|, |) vanishes at J' = 
irrespectively of N. Since the value of At(±|, |)/J' at J' = is almost 
independent of N, we can readily estimate this value in the limit of A^— ^oo 
to be 1.7 ± 0.1. Due to the finite-size effect, on the other hand, the value 
J'^Q^ of J' at which As(±|, |) vanishes is negative. This tends to as 
J'^Q^ ~ exp (—N/^') in the limit of A"— ^oo, where ^' is a constant. We have 
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estimated the infinite- value of As(±|, \)/ J' at J' = by extrapolating 
the finite- values of dAs(±|, |)/dJ' at J' = J'^q-^ to A^^oo. The result is 
1.1 d= 0.1. Thus, the value of the ratio at J' = obtained by the present 
numerical analysis is almost equal to |, which is in good agreement with 
the analytical result discussed above. When \J'\/J is small, the numerical 
calculation for the case of finite d also gives a satisfactorily good agreement 
with the analytical results. Details of the numerical calculation will be 
published in the near future. 

We summarize the results of the present study. 1) The analytical 
expressions for the energies of the low-lying states of H in the case of 
|J'|/J<Cl has been obtained by means of the variational method. 2) We 
have shown that 7i^^^ is equivalent to 7i in this case, as far as the energies 
of the low-lying states are concerned. 3) We have given the clear explana- 
tion for origin of the anisotropy of the exchange interaction in TY^^*^ [see 
eqs. (19) and (20)]. 4) The dependence of the energies on J' obtained by 
the variational method is in semi-quantitatively good agreement with that 
obtained by the numerical diagonalization. 
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Figure Captions 

Fig. 1. Plots versus doiAs{±^, \)/J' (dotted line), At(±|, f)/J' (dashed 
line), and At(±|, |)/J' (dot-dashed line) obtained for A = A' = 1 by 
the analytical method. 

Fig. 2. Schematic energy [measured from £^o,o+£t(=l=|5 |)] versus J' / J di- 
agram obtained for A = A' = 1 by the analytical method; (a) for (i>0, 
(b) for (i = 0, and (c) for d<Q. The dotted line shows As(±^, ^)/J, 
the dashed line shows At(±|, |)/J, and the dot-dashed line shows 
At(±^, |)/J. Note that in (b), At(±|, |)/J = At(±§, |)/J, which 
we show by the dashed line. The full line shows the bottom of the en- 
ergy continuum on the assumption that the value of the Haldane gap 
is equal to | J which is obtained by the analytical method for d = {)F) 

Fig. 3. Energy [measured from from £^o,o+£t(i^, |)] versus J' / J diagram 
obtained for A/" = 15 in the isotropic case of A = A' = 1 and <i=0 by the 
numerical diagonalization. The dotted line shows As(±|, \)/J and 
the dashed line shows At(±|, |)/J = At(±|, |)/J. 
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